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Abstract 

We study the ground state eigenvalues of Baxter's Q-operator for the eight-vertex model 
in a special case when it describes the off-critical deformation of the A = — ^ six- vertex 
model. We show that these eigenvalues satisfy a non-stationary Schrodinger equation with 
the time-dependent potential given by the Weierstrass elliptic p-function where the modular 
parameter r plays the role of (imaginary) time. In the scaling limit the equation transforms 
into a "non-stationary Mathieu equation" for the vacuum eigenvalues of the Q-operators in 
the finite- volume massive sine-Gordon model at the super-symmetric point, which is closely 
related to the theory of dilute polymers on a cylinder and the Painleve III equation. 
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1 Introduction and Summary 



The Q-operators introduced by Baxter in his pioneering paper on the eight-vertex model 
continue to reveal their exceptional properties in the theory of integrable quantum systems. 
These operators play a central role in the remarkable connection of Conformal Field Theory 
(CFT) with the spectral theory of the Schrodinger equation 2 discovered a few years ago jS]. 
As shown in [3115, the vacuum eigenvalues of the Q-operators j5| in CFT, with the central charge 
c < 1, can be identified with the spectral determinants of certain one-dimensional stationary 
Schrodinger equations. Some further developments and applications of this connection can found 
in the recent review article [0]. 

Apart from a few exceptions, the vacuum eigenvalues of the Q-operators (considered as func- 
tions of the spectral parameter) do not generally satisfy any ordinary second order differential 
equation themselves. One such exception is the case of the c = CFT where the these eigenval- 
ues for particular Virasoro vacuum states are known to satisfy the Bessel differential equation 
0. Remarkably, a similar property holds for the lattice counterparts of these eigenvalues in the 
A = — ^ six- vertex model (11111 for the chain of an odd number of sites. In this paper we explain 
how this property is generalized for the corresponding cases of the lattice eight-vertex model 
and the massive finite volume sine-Gordon model with N = 2 supersymmetry. 

We study the eight-vertex model on a periodic chain of an odd length, = 2n + 1, n = 
0, 1, 2, . . . oo. The eigenvalues of the transfer matrix of the model, T{u), satisfy Baxter's famous 
TQ-equation 

T{u) Q{u) = ^{u - T]) Q{u + 2r?) + + r/) Q{u - 2r/), (1) 
where u is the spectral parameter, 

0(n)=<(u|g), q = e'^\ Imr > 0, (2) 

and q) is the standard theta-function with the periods tt and vrr (we follow the notation 

of jni)- Here we consider a special case ry = vr/S, where the ground state eigenvalue is known 
jlll 1^1 to have a very simple form for all (odd) N 

T{u) = <P{u), r, = ^. (3) 

The equation (0) with this eigenvalue, T(u), has two different solutions El, Q±{u) = 
Q±iu,q,n), which are entire functions of the variable u and obey the following periodicity 
conditions [T1[T3]^ 

g±(u + 7r) = ±(-l)"Q±(u), Q±(n + vrr) =(7-^/2 e-*^" Q^(u), Q±{-u) = Q±{u). (4) 

The above requirements uniquely determine Q±{u) to within a common u-independent normal- 
ization factor. For further references it is convenient rewrite the functional equation for Q±{u) 
in the form 

27r 27r Att Att 

(^(n)Q±(n) + (j){u + —)Q±{u + y ) + + y )Q±(n + y) = 0- (5) 



We show that the functions 



*±(n) ^ ^±{u, q, n) = '^; Q±{u, q, n), (6) 



^The factor (-1)" in (gl and lUll reflects our convention for labeling the eigenvalues for different n, which will 
be important in Sec. 3 
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which are meromorphic functions of the variable u for any fixed values of q and n, satisfy the 
non-stationary Schrodinger equation 

d ( cP' ^ 

6g— ^'(n,g,n) = | - ^ + 9n(n + 1) ^(S'u | g^) + c{(i,n)y^[u,q,n), (7) 

where the modular parameter r plays the role of (imaginary) time and the time-dependent 
potential is defined through the elliptic Weierstrass p-function ^Oj (our function p(u | e*'^'^) 
has the periods vr and Tre). The constant c(g, n) appearing in (O is totally controlled by the 
normalization of Q±{u) and can be explicitly determined once this normalization is fixed (see 
Sec. 3 below). 

Equation is obviously related to the Lame differential equation and could be natu- 
rally called the "non-stationary Lame equation". To our knowledge this equation^ (in fact, 
a more general equation) first explicitly appeared in [T^ as a particular case of the Knizhnik- 
Zamolodchikov-Bernard equation |151 116j for the one-point correlation function in the sl(2)- 
WZW-model on the torus. Here, we will not explore this and other interesting connections 
of the equation ((TJ leaving that for the future. 

It is fairly trivial to show that the partial differential equation ((T)) for the meromorphic 
functions ^'±(u) is equivalent to the Baxter equation Q. Indeed, every solution of Q, with the 
required analytic properties in the variable u, implied by @ and satisfies eq.©. However, 
the very existence of exactly two solutions with these properties is by no means trivial and 
reflects some rather special features of eq.((7J discussed below. 

Equation Q has been discovered by virtue of a remarkable polynomial property of the eigen- 
values Q±{u). Using a combination of analytical and numerical techniques we have explicitly 
solved the equation @ for all values of n < 10. We have found that properly normalized 
eigenvalues Q±{u) could always be written as 

g±(n,g,n) = P^l \ g'^^M^ I , 7 = " 

where P2^^i{a, are homogeneous polynomials of the degree 2n + 1 in the variables a 
and P, with coefficients being polynomials in the variable 7 with integer coefficients. Then 
we considered a class of linear second order partial differential equations in two variables u 
and q where the coefficients of the second order derivatives are independent of n and all other 
coefficients are at most second degree polynomials in n. Eq.Q was then found as the only 
equation in this class satisfied by © with all explicitly calculated polynomials © with n < 10. 
It turned out that this equation uniquely defines two and only two such polynomials (jSJ for 
every value n = 0, 1,2, . . . ,00. It would be interesting to clarify the combinatorial nature of 
these polynomials, given that the related A = — | six- vertex model is connected to various 
important enumeration problems (18if ilQ]. 
In the scaling limit 

n^oo, q^O, t = 8ng^/^ = fixed, (9) 

the functions © essentially reduce to the ground state eigenvalues Q±{0) = Q±{9,t) of the Q- 
operators of the restricted massive sine-Gordon model (at the so-called, super-symmetric point) 

■^We are indebted to Prof. I.M.Krichever for informing us about the work |14|. 



^i(§|g^/^) 

^2(f kV2) 



(8) 



3 



on a cylinder of the spatial circumference R, where t = MR and M is the soliton mass. The 
equations (0) and (@J) become 

Q±{0) = Q±{e + 2TTi) + Q±{e-2TTi), (10) 

Q±{d + 37ri) = e^'^ Q±{e), Q+{0) = Q^{-e), (11) 

where the variable 9 is defined as u = ttt/2 — i9/3. With a suitable t-dependent normalization 
of Q±{0) equation Q could be brought to a particularly simple form 

Q±{e,t) = {^-^t^cosh2e-i)] Q±{e,t). (12) 

With the same normalization the asymptotic behavior of Q±{9) at large 6 is given by 

log Q± (^) = -\te' + log V±{t) + 2{dt log V±{t)-t/8) e-' 

+2(52logP±(t) -aaogP±(t)/t) e-2e + 0(e-3^), ^^+00, (13) 

where 'D±{t) are the Fredholm determinants which previously appeared in connection with the 
calculation of the "super symmetric index" and the problem of dilute polymers on a cylinder 
[1011211 muni- Note, in particular, that the quantity 

Fit) = 1 U(t). U(t) = log , (14) 

describes the free energy of a single incontractible polymer loop and satisfies the Painleve III 
equation [201 

In this connection, it is useful to mention the other celebrated appearances of Painleve transcen- 
dents in the theory of the two-dimensional Ising model and in the problem of isomonodromic 
transformations of the second order differential equations |25j . 

We did not attempt to make this paper self-contained. Detailed results are presented in |26j . 
The main reference for the eight-vertex model and its commuting T- and Q-matrices is Baxter's 
original paper Sec. 2 is meant to be read in conjunction with this paper. The definitions of 
the T- and Q-operators in continuous quantum field theory are given in |3 [23| . The connection 
of the Q-operators with the problem of dilute polymers is explained in |7]. 

2 The eight-vertex model and TQ-relation 

We consider the eight-vertex model on the A^-column square lattice with periodic boundary 
conditions and assume that is an odd integer = 2n + 1. Following we parameterize the 
Boltzmann weights a, b, c, d of the eight-vertex model as^ 

a = p ■di{2rj \ q^) ■^^{u -r]\q'^) t?i(n + 7/ 1 g^), 

b = p ■di{2r] I q^) 'di{u — 'q \ q^) t?4(n + rj \ q^), 

c = p ■di{2rj I q^) T?4(n — 'q\q^) i?4(n + r] \ q^), (16) 

d = p i)i{2r] I g^) 'di{u — r] \ q^) ■di{u + rj \ q^), 

^We use the notation of (lUI for theta-functions i?fc(it|g), k = 1,2,3,4, of the periods tt and ttt, q = e"^^, 
Imr > 0. The theta-functions H(t;), Q{v) of the nome qs used in ^ are given by 

qB = q\ H(«)=^i(^| g^), 0(„)=^4(^| 
where Kb is the complete elliptic integral of the first kind with the nome qs- 
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and fix the normalization factor p as 



p = 2 t?2(0|(z)-^ i94(0k')^^ (17) 
It is convenient to introduce new variables 7 and x, where 7 is defined by (cf. PP), 



{a — h + c — d){a — h — c + d) 

7 ^ 



{a + h + c + d){a + h-c-d) '92{r]\q^/'^) 
and x satisfies the quadratic equation 

7 \2 16 (a- 6)2 cd 



(18) 



(19) 



y/xJ {c + d)"^ {a + h + c + d)[a + h- c- d)' 

where we choose the following root 

2 

x = 7^, Un)=U^W/% Mn)=U^\q"'). (20) 
<(w) 2 2 

The row-to-row transfer-matrices T(u) and the Q(ti)-matrices of the model form a com- 
mutative family; their eigenvalues satisfy the functional relation In ^ Baxter explicitly 
constructed the matrix Q(ti) provided that 77 = ^ for integer m and L. We only consider a 
special case when the weights (|16() satisfy the constraint IP 

(a^ + ab) (62 + a6) = (c^ + a6) (d^ + ah), (21) 

which is equivalent to the condition = 7r/3. In jHl EHI it was conjectured that the largest 
eigenvalue of the transfer matrix (corresponding to the double-degenerate ground state of the 
model) has the simple form, 

T{u) = {a + h)^ = (t){u). (22) 

Here we study the corresponding eigenvalues Q±{u) of the Q-matrix. As noted in the 
method used in ^ for the construction of the Q-matrix cannot be executed in its full strength for 
7] = vr/S, since some axillary Q-matrix, CIr{u), in is not invertible in the full 2''^-dimensional 
space of states of the model. The numerical results presented in Table 1 of [IH] make it natural to 
suggest that the rank of Q,r{u) in this case is given by the A^-th Lucas number ((1 + \/5)/2)^ + 
((1 — -v/5)/2)^ which coincides with the dimension of the space of states of the A^-site hard 
hexagon model jSHj^- This indicates that for 77 = 7r/3 the construction of 1 only provides a 
"restricted" Q(^)(ti)-matrix which acts in some "RSOS-projected" subspace of the full space 
of states of the eight-vertex model. This interesting phenomenon certainly deserves special 
investigations in its own right; we have verified for several small values of A^ that the ground 
state eigenvectors corresponding to (|22|) belongs to this RSOS-projected subspace and that the 
eigenvalues of Cl^'^\u) satisfy Below we will assume that this is true for all (odd) A^. 

After this brief review, let us now describe our main results. For r/ = vr/S the variable 
7 = 7(g), defined by H18() . depends on q only, while the variable x = x{u,q), defined by (fT^. 
depends on u and q. Below it will be more convenient to use the combinations 

Qi{u) = {Q+{u) + Q-{u))/2, Q2{u) = {Q+{u) - Q^{u))/2. (23) 



^S.M.Sergeev noted jl that only a minor modification of the arguments of ^ leads to the matrix with 
the rank equal to 2^ for even iV and to (2^^ - 2) for odd A^. 
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which are simply related by the periodicity relation 

Qi,2{u + n) = {-irQ2,i{u). (24) 

Bearing in mind this simple relation we will only quote results for Qi{u), writing it as Q'f'\u) 
to indicate the n-dependence. We have found that all the eigenvalues Q^^\u) can be written as 

Qf\u)=M{q,n)l^{u)dl''{u)Vn{x,z), z = (25) 

where J\f{q,n) is an arbitrary normalization factor and Vnix,z) are polynomials in x,z of the 
degree n in a;, 

n 

in) I 

k=0 



Vn{x,z) = Y,rt"\z)xK (26) 



while r^"\z), i = 0, ...,n, are polynomials in z with integer coefficients. The normalization 
of Vn{x,z) is fixed by the requirement ri"^(0) = 1. The polynomials Vn{x,z) are uniquely 
determined by the following partial differential equation in the variables x and z, 

\^A{x,z)dl + Bn{x,z)d^ + Cn{x,z) + T{x, z) 8,} Vn{x, z) = 0, (27) 

where 

A{x,z) = 2x{l + x-3xz + x^z){x + Az-6xz-3xz^ + Ax^z^), (28) 

Bn{x,z) = A{l + x-3xz + x^z){x + ?>z-7xz + ?>x^z^)+ (29) 

+2nx(l - 14z + 21^2 - 8x^z^ + 3x^2(32^ + 6z - 1) - x(l - 9z + 23z^ + 9z^)) 
Cn{x,z) = n[z(9z - 5) + x2z(3z2 + llz - 2) + x(9z^ - 38^2 + 192 - 2) - 

-4x^z^ + nz{l -9z- x{9z^ - 36z + 3) + x^iSz"^ - 31z + 4) + 8x^2^)] (30) 

T{x,z) = -2z(l - z)(l - 9z)(l + x - 3xz + x^z). (31) 

The first few polynomials Vnix, z) read 

Vq{x,z) = 1, Vi{x,z) = x + 3, V2{x,z)=x^{l + z) + bx{l + 3z) + l{), (32) 

7^3(2;, z) = x^(l + 3z + Az^) + 7x^(1 + 5z + 18^^) + 7x(3 + 19z + ISz^) + 35 + 21z. (33) 

In the next Section, we will prove that the eigenvalues Qi""* i'^) given by (|25j) (as well as all related 
eigenvalues Q^2\''A ^'^^ Q± '*('") given by (|23|) and ((SU) automatically satisfy the functional 
relation ©, merely as a consequence of the defining property H27() of Vn{x.,z). Of course, for 
small values of n this functional relation can be checked directly. For example, it is not very 
difficult to check it for 

Qf\u) ~ ^3(n), Q^i\u) ~ ^3(n) [i^l{u) + 3 9l{u)\ , (34) 

by employing various identities for elliptic functions, however for n>2 this does not appear to 
be practical. 

The polynomials Vn{x,z) can be effectively calculated with the following procedure. It is 
easy to see that (|27|) leads to descending recurrence relations for coefficients in (|26j) in the sense 
that each coefficient r'^j^\z) with k < n can be recursively calculated in terms of rm\z), with 
m = k + 1, . . . ,n and, therefore, can be eventually expressed through the coefficient Sn{z) = 
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rli'\z) of the leading power of x. These leading coefficients, Sn{z), n = 0, 1, 2 . . . oo, are uniquely 
determined by the following recurrence relation 

2z(z - l){9z - l)2[logs„(z)];' + 2(32 - lf{9z - l)[log s„(z)];+ (35) 
+8(2n + i)2 ^n+i(^K-i(^) _ ^4^3^ ^ ^^^3^ ^ 2) + {9z - l)n(5n + 3)] = 0, 

with the initial condition So{z) = si{z) = 1. In particular, for z = 1/9 (corresponding to q = 0) 
this gives 

2"-in!(3n + 2)! 
3"((2n + 1)!)2 

Currently, the polynomials Vn{x,z) have been calculated explicitly for n < 50; using the above 
procedure they can be easily calculated for higher values of n. 



50(9) = 1, Sn+l{g) = onU0^^^m2 ^n{g). (36) 



3 Non- stationary Lame equation 

Evidently, the algebraic form of the partial differential equation given by ()27|) is not very illumi- 
nating (even though it is quite useful for the analysis of polynomial solutions). Fortunately this 
equation has a much more elegant equivalent form given by the non-stationary Lame equation 
(|7|) discussed in the Introduction. The details of transformations between these two forms will 
be presented elsewhere |26j. Noting the frustrating expressions for the coefficients in H27|) it is 
not surprising that these transformations turn out to be rather tedious. They involve many 
elliptic function identities, in particular, the algebraic properties of the ring of theta-constants 
and their (^-derivatives |32| happened to be extremely useful. 
Let us choose the normalization factor Af{q, n) in as 

n(n + l) 

(-1)" 2(3"'+5"+i)/2 ) Ar(g, n) = 7"" ^[{0 I g)-"-V3 |(^2 _ ^4(0 I 2 ^3^^ 

and define the functions \I'±(u, g, n) as in © where Q±{u) are given by (|23|) . (|2H) and (|25|) . 
Then the equation (|27|1 takes the form ((T)) where the constant term c{q, n) is given by 

c{q,n) = -2,n{n + l)^j^^. (38) 

By construction the functions ^±{u,q,n) are meromorphic functions of the variable u, which 
obey periodicity relations 

^(u + 27r) = -^{u), -^{u + 2-kt) = q'^ e"^*" ^(u), -^{-u) = (-1)'^+i^'(m) (39) 

and have {n + l)-th order zeros a,t u = kn + rmrT, where k,m £ i.e., 

^'(e + Arvr + mvrr) = 0(e"+^), e ^ 0, k,m£Z. (40) 

Let us show that for n > 1, equation 0, restricted to a class of functions ^'(m) with such 
analytic properties, is equivalent to the functional relation Q. Any such solution of Q, ^'(m), 
could have either an (n + l)-th order zero or an n-th order pole in the variable u, at all points 
u = {3k ± l)'/r/3 + mvrr, with k,m G and these are the only points where ^{u) could have 
poles. Consider the function 

2tt Ah 

^{u) = ^{u) + ^{u + —) + ^{u + —), (41) 
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which also satisfies ((7J along with ^'(n). When u = the second and third terms in l|41j) may 
have n-th order poles, however, they must cancel each other due to the last relation in Thus 
could have at most (n— l)-th order pole at n = (the first term in (|41|) obviously vanishes 
due to (jini)). As noted above such pole is forbidden by eq.© and thus ^{u) has the (n + l)-th 
order zero at n = 0. Similarly, one concludes that vanishes at all points u = kir/S + mvrr, 
with k,ni £ 7j, and, therefore, has at least 12(n + 1) zeroes in the periodicity parallelogram (of 
the periods 2tt and 2ttt) and no poles at all. However, for n > 1 this contradicts (|39|) . unless 
$(?/) = 0, which is equivalent to ©. The special case n = is considered in 

Now consider various limiting forms of the equation ((T)). When q ^ 0, with u and n kept 
fixed, the functions ^±{u) reduce to their analogs for the 6- vertex model 

^±{u,q,n)=q'2^'^±+i^^^l''\u,n){l + 0{q)), q ^ 0, (42) 

where d-t = (1 =F 6)/36 and eq.© becomes 

which is simply related to eq.(13) of 9 . Taking now n ^ oo and s ~ ilogn one recovers the 
eigenvalues of the Q-operators of ref. [5^ corresponding to p = ±1/6 vacuum states in the c = 
CFT 

q{cft) ^ ^ J ^.^ ^(6^) . i iog(8„/i) _ iQ\ ^ (44) 

n— >oo 

which are known to satisfy the Bessel differential equation 

{-dl + de + ^t^ e^' + d±]Q^^''^\e + logt) = 0. (45) 
In the most interesting scaling limit Q, the limiting functions 

Q±(e, t) = t^/^ e*'/i6 ^-1/4 ^^(7rr/2 - e^"", n) (46) 



-i:log(8n/t) 



coincide with the eigenvalues |7j of the Q-operators for special twisted vacuum states in the 
massive sine-Gordon at the supersymmetric point j22j (where the ground state energy of the 
model vanishes identically due to the supersymmetry) . These eigenvalues satisfy the "non- 
stationary Mathieu equation" ()12() . Let us show that this equation completely determines the 
asymptotic expansion of these eigenvalues at large 0, 



log Q±{e,t) = --/ + Y,S±it) e'^', e ^ +00. (47) 

fc=0 



Consider the integral operator fC{t) with the kernel 



I -u{e)-u(e') ^ 
mO,e') = — , n(0) = -coshe, (48) 



which satisfies the following identity 



tdt+Me-Me' lC{t\e,e') = —te-<^^-''^^'^+^\ (49) 

47r 
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where Mg denotes the differential operator in the RHS of Using this identity one can show 
that Q±{9) satisfy the hnear integral equation discovered in [7j 

/CO 
ic{t\e,e') Q±{e\t)de\ (so) 
-oo 

provided that the functions T)±{t) coincide with the Fredholm determinants |22j 

V±{t) =C±dQl{l± K:{t)), (51) 

where C± are numerical constants. Comparing H47|) with (jSUf) one concludes that {t) = 
logP±(i). Then eq. (|12|) allows to express all higher coefficients in (|47|) through the power 
{t) and its derivatives. The few first coefficients are shown in (|13() . Finally note, that in the 
limit t ^ 0, 6* ~ - log t, 

Q±{e,t) = t<^±Qf^^\e + \ogt){i + o{t^/^)), t^o, e^-iogt, (52) 

and the equation (fT^ reduces to as it, of course, should^. 

We expect that our results can be readily extended to elliptic generalizations of the other 
special lattice models |33| . closely related to the rj = tt/S six- vertex model. However, it would 
be more important to understand whether similar considerations could be applied for the eight- 
vertex with an arbitrary value of r] and whether the general scheme of correspondence between 
the c < 1 CFT and ordinary differential equations, developed in [HlllJISJI, can be extended to 
the massive field theory (sine-Gordon model) with the use of partial differential equations. The 
investigation of these questions is in progress. 
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